9 Search Results

Improved confinement in negative triangularity (NT) experiments is attributed to reduced fluxes driven by micro-turbulence. The physical mechanism of why thermal confinement improves in NT relative to PT is unknown. This study employs gyrokinetic flux tube simulations using the GENE code with local Miller equilibrium to elucidate the physical mechanisms behind the beneficial effects of NT flux surface shapes. The focus is on collisionless ion temperature gradient (ITG) driven turbulence with adiabatic electrons. The kinetic profiles are held fixed across a scan of triangularity values, thus enabling comparisons on a level playing field. The reduced linear growth rates for NTmore » is shown to be due to a reduced eigenmode averaged magnetic drift frequency and a wider, stronger negative local magnetic shear region about the outboard mid-plane. The nonlinear heat flux is lower for NT than that for PT, due to reduced radial correlation length and increased correlation time (τ_c) of fluctuations. These, in turn, are due to a comparatively higher level of self-generated zero-frequency E × B zonal shearing rate ω_E in NT as compared to PT. Though the linear zonal potential residual is lower for NT, the nonlinearly generated E × B zonal shearing rate is higher for NT than for PT. This outcome is linked to the distinctive features of the radial wavenumber spectra of the zonal potential and the zonal shearing rate. The dimensionless parameter ω_Eτ_c is suggested as a figure of merit. This is higher for NT than for PT. Thus, the reduced heat diffusivity for NT is linked to increased ω_Eτ_c. Self-generated temperature corrugations (i.e. zonal temperature gradients) are much weaker than the background mean temperature gradient. Nevertheless, temperature corrugations are more pronounced in NT than in PT.« less

Abstract This paper presents a comparative study of the poloidal distribution of the mean E  ×  B shearing rate for positive triangularity (PT) and negative triangularity (NT) tokamaks. The effects of flux surface up–down asymmetry due to asymmetric upper and lower triangularities are also considered. Both direct eddy straining and the effects on Shafranov shift feedback loops are examined. Shafranov shift increases the shearing rate at all poloidal angles for all triangularities, due to flux surface compression. The maximum shearing rate bifurcates at a critical triangularity ${\delta }_{\text{crit}}(\lesssim 0)$more » . Thus, the shearing rate is maximal off the outboard midplane for NT, while it is maximal on the outboard midplane for PT. For up–down asymmetric triangularity, the usual up–down symmetry of the shearing rate is broken. The shearing rate at the outboard midplane is lower for NT than for PT, suggesting that the shearing efficiency in NT is reduced. Implications for turbulence stabilisation and confinement improvement in high- ${\beta }_{\mathrm{p}}$ NT and internal transport barrier discharges are discussed.« less

Here, this paper presents a comparative study of zonal flow screening in positive and negative triangularity tokamaks. Neoclassical screening sets the strength of zonal flow shear for fixed drive and damping. Orbit calculations show that the radial excursions of trapped particles are reduced in negative triangularity tokamaks, as compared to positive triangularity. Yet surprisingly, the neoclassical dielectric susceptibility actually increases with decreasing triangularity, such that the susceptibility is higher for negative triangularity than for positive triangularity. This is because the reduction in neoclassical polarization by shrinking the banana width is offset by the increase in neoclassical polarization by the enhancementmore » of trapped fraction for negative triangularity. As a result, the zonal flow screening length is actually enhanced for negative triangularity, as compared to positive triangularity. Hence, the zonal flow residual is smaller for negative triangularity than for positive triangularity. Results from gyrokinetic simulations support the analytic calculations.« less

Electric field profile structure—especially its shear—is a natural order parameter for the edge plasma, and characterizes confinement regimes ranging from the H-mode (Wagner et al. 1982 Phys. Rev. Lett.49, 1408–1412 (doi:10.1103/PhysRevLett.49.1408)) to the density limit (DL) (Greenwald et al. 1988 Nucl. Fusion28, 2199–2207 (doi:10.1088/0029-5515/28/12/009)). The theoretical developments and lessons learned during 40 years of H-mode studies (Connor & Wilson 1999 Plasma Phys. Control. Fusion42, R1–R74 (doi:10.1088/0741-3335/42/1/201); Wagner 2007 Plasma Phys. Control. Fusion49, B1–B33 (doi:10.1088/0741-3335/49/12b/s01)) are applied to the shear layer collapse paradigm (Hong et al. 2017 Nucl. Fusion58, 016041 (doi:10.1088/1741-4326/aa9626)) for the onset of DL phenomena. Results from recent experimentsmore » on edge shear layers and DL phenomenology are summarized and discussed in the light of L → H transition physics. The theory of shear layer collapse is then developed. In this work, we demonstrate that shear layer physics captures both the well known current (Greenwald) scaling of the DL (Greenwald 2002 Plasma Phys. Control. Fusion44, R27–R53 (doi:10.1088/0741-3335/44/8/201); Greenwald et al. 2014 Phys. Plasmas21, 110501 (doi:10.1063/1.4901920)), as well as the emerging power scaling (Zanca, Sattin, JET Contributors 2019 Nucl. Fusion59, 126011 (doi:10.1088/1741-4326/ab3b31)). The derivation of the power scaling theory exploits an existing model, originally developed for the L → H transition (Diamond, Liang, Carreras, Terry 1994 Phys. Rev. Lett.72, 2565–2568 (doi:10.1103/PhysRevLett.72.2565); Kim & Diamond 2003 Phys. Rev. Lett.90, 185006 (doi:10.1103/PhysRevLett.90.185006)). We describe the enhanced particle transport events that occur following shear layer collapse. Open problems and future directions are discussed.« less

Edge shear layer collapse causes edge cooling and aggravates radiative effects. This paper details on the microscopic dynamics of the emergence of power (Q) scaling of density limit (DL) from the shear layer collapse transport bifurcation scenario. The analysis is based on a novel 4-field model, which evolves turbulence energy, zonal flow energy, temperature gradient and density, including the neoclassical screening of zonal flow response. Bifurcation analysis yields power scaling of critical density for shear layer collapse as $$n_{crit}\sim Q^{1/3}$$. The favorable Q scaling of the DL emerges from the fact that the shear layer strength increases with Q, thusmore » preventing shear layer collapse. This in turn reduces particle transport and improves particle confinement. RMP induced ambient stochastic fields degrade the shear layer by inducing decoherence in the Reynolds stress. As a result the particle transport increases and particle confinement degrades. This leads to the emergence of unfavorable stochastic field intensity ($$b_{st}^{2}$$) scaling of the critical density as $$n_{crit}\sim(1+b_{st}^{2})^{-5/3}$$. All fields, including zonal flow shear, exhibit hysteresis when the power (Q) is ramped cyclically across the bifurcation point. The hysteresis is due to dynamical delay in bifurcation on account of critical slowing down. Furthermore, the dynamical hysteresis here is fundamentally different from the hysteresis associated with the existence of bi-stable states.« less

This paper details the theory of edge shear layer collapse as the density approaches the Greenwald density limit. It significantly extends earlier work, which was restricted in applicability. The zonal shear flow screening length is calculated for banana, plateau and Pfirsch–Schluter regimes. Poloidal field scaling persists in the plateau regime. Neoclassical screening and drift wave–zonal flow dynamics are combined in a theory, which is then reduced to a predator–prey model. Zonal noise, due to incoherent mode coupling, is retained. The threshold condition for edge shear layer collapse is computed, and linked to a critical value of the dimensionless parameter $${\rhomore » }_{\mathrm{s}}/\sqrt{{\rho }_{\mathrm{s}\mathrm{c}}{L}_{n}}$$. Here $${\rho }_{s}$$ is the ion sound radius, $${\rho }_{sc}$$ is the zonal flow screening length and $${L}_{n}$$ is the equilibrium density scale length. The limiting initial edge density for shear layer collapse is derived and shown to scale favorably with the plasma current. Here, the results are discussed in light of the density limit and Ohmic phenomenology.« less

A unified theory of zonal flow shears and density corrugations in drift wave turbulence is presented. Polarization and density advection beat excitation are studied in combination with modulational response. Noise is driven by two-time flux correlation. While the effective zonal flow eddy viscosity can go negative, the zonal diffusivity is positive definite. There is no inverse cascade of density corrugation. The connection between avalanches and corrugations is discussed. Here, the zonal cross-correlation is identified and calculated. Conditions for alignment of zonal shears and corrugation gradients are determined, and the implications for staircase structure are discussed. We show that the synergymore » of beat noise and modulational effects is stronger than either alone. Strong zonal flows can be excited well below the modulational instability threshold. In the context of L–H transition, zonal noise quenches turbulence overshoot by eliminating the threshold for zonal flow excitation. The power threshold for L–H transition is lowered.« less

Tangled magnetic fields, often coexisting with an ordered mean field, have a major impact on turbulence and momentum transport in many plasmas, including those found in the solar tachocline and magnetic confinement devices. We present a novel mean field theory of potential vorticity mixing in β-plane magnetohydrodynamic (MHD) and drift wave turbulence. Our results show that mean square stochastic fields strongly reduce Reynolds stress coherence. This decoherence of potential vorticity flux due to stochastic field scattering leads to suppression of momentum transport and zonal flow formation. A simple calculation suggests that the breaking of the shear-eddy tilting feedback loop bymore » stochastic fields is the key underlying physics mechanism. Furthermore, a dimensionless parameter that quantifies the increment in power threshold is identified and used to assess the impact of stochastic field on the L-H transition. We discuss a model of stochastic fields as a resisto-elastic network.« less

Few, if any, of the many papers on turbulence spreading address the key question of how turbulence spreading actually affects the profile structure. In this work, we are using a reduced model to answer that question. Turbulence spreading is most relevant near regions where the profiles support a strong intensity gradient ∇I. One such case is at the edge of an L mode ischarge, near a source of turbulence [i.e., either a localized source of edge turbulence or an influx of turbulence from the scrape-off layer (SOL)]. Another is in “No Man’s Land” (NML), which connects the pedestal to themore » stiff core in H mode. In the case of L mode, without an edge intensity source, the turbulence intensity profile is nearly flat and spreading has a weak effect. An edge localized source increases the edge ∇I, which then drives inward spreading. Invasion of turbulence from the SOL to the edge softens the edge pressure gradient. In H mode, the strong shear suppression of pedestal turbulence necessarily forces a sharp ∇I in NML. This sharp∇I drives a significant flux of turbulence from the core to the pedestal, where it is ultimately dissipated by shearing. Counter-intuitively, the results indicate that spreading actually increases the pedestal height and width and hence the energy content in H mode. This suggests that models of the pedestal structure should include NML turbulence spreading effects. The relation of avalanches to spreading is studied. Spreading weakly affects the avalanche distribution, but the spatiotemporal correlation of intensity increases with spreading.« less